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ABSTRACT 

In the last years, more and more interest has been devoted to analytical solutions, 
including inflow and outflow, to study the metallicity enrichment in galaxies. In this 
framework, we assume a star formation rate which follows a linear Schmidt law, and 
we present new analytical solutions for the evolution of the metallicity (Z) in galaxies. 
In particular, we take into account environmental effects including primordial and 
enriched gas infall, outflow, different star formation efficiencies and galactic fountains. 
The enriched infall is included to take into account galaxy-galaxy interactions. Our 
main results can be summarized as: i) when a linear Schmidt law of star formation 
is assumed, the resulting time evolution of the metallicity Z is the same either for a 
closed-box model or for an outflow model, ii) The mass-metallicity relation for galaxies 
which suffer a chemically enriched infall, originating from another evolved galaxy with 
no pre-enriched gas, is shifted down in parallel at lower Z values, if compared the 
closed box model, hi) When a galaxy suffers at the same time a primordial infall 
and a chemically enriched one, the primordial infall always dominates the chemical 
evolution, iv) We present new solutions for the metallicity evolution in a galaxy which 
suffers galactic fountains and an enriched infall from another galaxy at the same 
time. The analytical solutions presented here can be very important to study the 
metallicity (oxygen), which is measured in high-redshift objects. These solutions can 
be very useful: a) in the context of cosmological semi-analytical models for galaxy 
formation and evolution, and b) for the study of compact groups of galaxies. 
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1 INTRODUCTION 

The galactic chemical evolution is the study of the trans¬ 
formation of gas into stars and the resulting evolution of 
the chemical composition of a galaxy. The so called “Sim¬ 
ple Model” remains an useful guide for understanding the 
chemical evolution of galaxies since the pioneering works of 
Schmidt (1963), Searle & Sargent (1972), Tinsley (1974), 
Pagel & Patchett (1975). In order to derive analytical so¬ 
lutions for the chemical evolution of galaxies one needs to 
make several hypotheses: the initial mass function (IMF) 
should be considered constant, the lifetime of stars should 
be neglected (instantaneous recycling approximation, IRA), 
and the complete mixing of chemical elements with the sur¬ 
rounding interstellar medium (ISM). These assumptions, in 
fact, allow us to have analytical expressions for the metallic¬ 
ity evolution of the galaxies in time and in terms of the gas 
fraction. One can find analytical solutions even in presence 
of gas flows (infall, outflow, galactic fountains) as shown 
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by Tinsley (1980), Clayton (1988), Lacey & Fall (1985), 
Matteucci & Chios! (1983), Edmunds (1990), Recchi et al. 
(2008), Spitoni et al. (2010), Peeples & Shankar (2011), Lilly 
et al. (2013), Pipino et al. (2014), Peng & Maiolino (2014a), 
Recchi & Kroupa (2014), Kudritzki et al. (2015). 

However, in all mentioned works the solutions are ob¬ 
tained only under specific assumptions regarding the in¬ 
fall/outflow rates. However, we note that analytical solutions 
are not able to give a complete description of the chemistry 
of a galaxy, since they fail in following the evolution of el¬ 
ements created on long time-scales, such as iron and nitro¬ 
gen. A satisfactory description of the iron evolution requires 
detailed numerical models relaxing IRA and including the 
chemical enrichment from Type la SNe (e.g. Matteucci & 
Greggio 1986) allow one to follow in detail the evolution of 
single elements. Matteucci & Greggio (1986) showed in de¬ 
tail the effect of the time-delay model, already suggested by 
Tinsley (1980) and Greggio & Renzini (1983), in particular 
the effect of a delayed Fe production by Type la SNe on 
abundance ratios involving a-elements (O, Mg, Si). Analyt¬ 
ical solution can be, on the other hand, adopted when study- 
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ing the evolution of oxygen, created on short timescales and 
tracing the evolution of the global metallicity, Z, of which 
oxygen is the main component. 

In this work we start by using the formalism described 
by Matteucci (2001), Recchi et al. (2008), and Spitoni et al. 
(2010) and we show, for the first time, the analytical solution 
of the evolution of the metallicity of a galaxy in presence of 
“environment” effects coupled with galactic fountains and 
primordial infall of gas. In this context for “environment” 
effect we mean the situation where a galaxy suffers, during 
its evolution, infall of enriched gas from another evolving 
galactic system, and this gas represents an enriched infall 
variable in time. 

The dynamics of interacting systems have been the 
subject of many papers concerning numerical simula¬ 
tions (Toomre & Toomre 1972; Barnes & Hernquist 1992; 
Berentzen et al. 2003) and spectrophotometric models (Lar¬ 
son & Tinsley 1978; Kennicutt 1990; Temporin et al. 2003a, 
2003b). Large amounts of ISM can be removed from the 
main disks of spiral galaxies by different processes: tides due 
to the gravitational force of a companion, ram pressure strip¬ 
ping during a near head-on collision between gas-rich galax¬ 
ies, ram pressure stripping by intracluster gas, and galactic 
winds driven by supernovae. Smith & Struck (2001) using 
CO signatures observed 11 extragalactic tails and bridges 
in nine interacting galaxy systems. Recently, Smith et al. 
(2010) using Galaxy Evolution Explorer (GALEX) ultravio¬ 
let telescope studied star formation morphology and stellar 
populations and in 42 nearby optically selected pre-merger 
interacting galaxy pairs. Tails and bridges structures are of¬ 
ten more prominent relative to the disks in UV images com¬ 
pared to optical maps. This effect is likely due to enhanced 
star formation in the tidal features compared to the disks 
rather than to reduced extinction. We also refer the reader 
to the review of Boselli & Gavazzi (2006) where a compre¬ 
hensive description of the environment effects on late-type 
galaxies in nearby clusters is presented. Among them we re¬ 
call the tidal interactions among galaxy pairs which act on 
gas, dust and stars, as well as on dark matter and is depend¬ 
ing on the gravitational bounding of the various components; 
tidal interaction between galaxies and the cluster potential 
well, and finally the so-called “galaxy harassment”: the evo¬ 
lution of cluster galaxies is governed by the combined effect 
of multiple high speed galaxy-galaxy close (~ 50 kpc) en¬ 
counters with the interaction with the potential well of the 
cluster as a whole. 

As stated by Davies et al. (2010), galaxy-galaxy inter¬ 
actions are particularly important and striking when the 
speed of the interaction is well matched to the velocities 
of the stars and gas. Therefore, small galaxy groups can po¬ 
tentially provide the environment for dramatic gravitational 
disturbances. Beyond the Local Group, the closest example 
of this is the environment around M81, the M81 group. In 
fact, in M81 group it is evident that various galaxies are 
connected by flows. Extended filamentary structures exter¬ 
nal to the disc of M81, are clearly seen in emission in all of 
the Herschel bands. These complex interactions cannot be 
described by the simple tools offered, for example, by Rec¬ 
chi et al. (2008) or similar papers in literature, and with this 
paper we provide a new set of analytical solutions to take 
into account this kind of interactions in the framework of 
analytical chemical evolution models. 


Moreover, an enriched infall of gas can be originated 
from the galaxy itself. In fact, in the galactic fountain mod¬ 
els (Shapiro & Field 1976, Houck & Bregman 1990), hot 
gas is ejected out of the Galactic disk by supernova (SN) 
explosions, and part of this gas falls back in the form of 
condensed neutral clouds which move at intermediate and 
high radial velocities. For example, in the galactic fountain 
model the ejected gas from SN events falls back ballistically 
(Bregman 1980). These models are able to explain the ver¬ 
tical motion of the cold and warm gas components observed 
in several spiral galaxies (e.g. Fraternali et al. 2004; Boosma 
et al. 2005). Following the analytical implementation of the 
galactic fountain presented by Recchi et al. (2008), we will 
consider such an effect in our new analytical solutions. 

In the context of cosmological semi-analytical models of 
galaxy formation and evolution each galaxy is treated as one 
unresolved object, using integrated properties to describe 
the mass of stars, cold gas, hot gas and the black hole. Since 
each component of the galaxy is represented by one number, 
the dynamics within the galaxy is not resolved, and one 
needs to assume with laws for star formation, cooling and 
feedback that are valid on average for the entire galaxy. In 
this context, our new analytical solutions shall be extremely 
useful, because they give simple recipes concerning the time 
evolution of the global metallicity of a galaxy in different 
situations including inflow and outflow. 

The paper is organized as follows: in Sect. 2 we present 
the main assumptions of the “analytical” chemical evolu¬ 
tion models, in Sect. 3 we describe our assumptions and the 
system of equation we need to solve, and in Sect 4 new an¬ 
alytical solutions are presented. Finally, our conclusions are 
summarized in Section 5. In Appendix A we draw the com¬ 
plete expressions of some new analytical solutions we pre¬ 
sented in this paper, and in Appendix B a list of variables 
and parameters used throughout the paper is reported. 

2 THE CLOSED BOX AND LEAKY-BOX 
PRESCRIPTIONS 

The main assumptions of the Simple Model (Tinsley 1980) 
are: 

(i) The IMF is constant in time. 

(ii) The gas is well mixed at any time (instantaneous mix¬ 
ing approximation). 

(iii) Stars ^ 1 M© die instantaneously; stars smaller than 
1 M© live forever (instantaneous recycling approximation or 
IRA). 

These simplifying assumptions allow us to calculate analyt¬ 
ically the chemical evolution of the galaxies once we have 
defined the fundamental quantities, such as the returned 
fraction: 

/ OO 

(m — Mii)(f)(m)dm, (1) 

(where <j)(m) is the IMF and Mn is the mass of the remnant) 
and the yield per stellar generation: 

1 f°° 

yz = YTTr J rnpz,m4’(ni)dm, (2) 
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(where pz,m is the fraction of newly produced and ejected 
metals by a star of mass m). 

Recently, Recchi & Kroupa (2015) applied the inte¬ 
grated galactic initial mass function (IGIMF) to the sim¬ 
ple model solution, and yz and R are not constant, but are 
functions of time,through the time dependence of the star 
formation rate (SFR) and the metallicity of the system. 

The well known solution of the so called “closed-box” 
where the system is one-zone and there are no inflows nor 
outflows with constant mass (gas plus stars) is: 

Z = yz\n{fi~^) (3) 

where y is the gas fraction Mgas/Mtot, with Mtot = M, -|- 
Mgas- It is also assumed that the Mgaaifi) = Mtot{0) and 
the initial metallicity of the system is zero. 

Analytical solutions of simple models of chemical evo¬ 
lution including infall or outflow are known since at least 30 
years (Pagel & Patchett 1975, Hartwick 1976, Clayton 1988, 
Twarog 1980, Edmunds 1990). Here, we follow the approach 
and the terminology of Matteucci (2001), where it was as¬ 
sumed for simplicity linear flows ( gas flows are proportional 
to the SFR). Therefore, the outflow rate W{t) is defined as: 

W{t) = X{l-R)m, ( 4 ) 

where i/>(t) is the SFR, and the infall rate A{t) is given by: 

A{t) = A(1 - R)V>(t). (5) 

Here A and A are two proportionality constants ^ 0. The 
first assumption is justified by the fact that, the larger the 
SFR is, the more intense are the energetic events associated 
with it (in particular supernova explosions and stellar winds) 
and therefore the larger is the chance of having a large-scale 
outflow (see e.g. Silk 2003). A proportionality between A(t) 
and '0(t) has been discussed by Recchi et al. (2008). They 
tested different prescriptions for the infall of gas showing 
that their results do not change substantially if a generic ex¬ 
ponential infall is assumed. This demonstrates also that the 
major source of error in the solutions of the simple models 
is the IRA assumption rather than the assumption of linear 
flows. 

In Lilly et al (2013) and Pipino et al. (2014), they con¬ 
sidered the simple chemical evolution model in the cosmo¬ 
logical context where the gas accreted is proportional to the 
dark matter growth. In particular, in Pipino et al. (2014) 
the infall parameter A(t) is time dependent and is defined 
as the ratio between the accretion rate from cosmological 
simulations and the SFR. In this paper we do not inves¬ 
tigate “cosmological” aspects but we study the effects of 
galaxy-galaxy interactions on the time evolution of galactic 
metallicity. 

The evolution of the metallicity of a system as a func¬ 
tion of /i in the case with only outflows, i.e. A(t) = 0 and 
W{t) ^ 0 is the following: 

(T^H(1 + A)m''-A]. (6) 

with the assumption that at t = 0, Z(0)=0, Mtotifl) = 
Mgas{Q). In the opposite case {A{t) / 0 and W{t) = 0), 
assuming for the infalling gas a primordial composition (e.g 
Za = 0), we obtain this solution (Matteucci 2001): 

( 7 ) 


The general solution presented by Matteucci (2001) for a 
system described by the simple model in the presence of 
infall of gas with a general metallicity Za and outflow is: 

^ = ^ZA + yz 1 1 ^ I ( 8 ) 


3 THE CHEMICAL EVOLUTION MODEL 

WITH ENVIRONMENT EFFECTS 

Using the formalism introduced by Matteucci (2001), and 
adopted later by Recchi et al. (2008) and Spitoni et al. 
(2010), for the hrst time we intend to study the effects of 
the environment on the chemical evolution of a galaxy by 
means of analytical solutions. The environment effects are 
mimicked by an infall of gas with a time dependent metal¬ 
licity originated by a nearby galaxy. 

In our model we analyze the evolution of the oxygen in 
a generic galaxy, and assume for the oxygen yield standard 
values of yo = 0.01 and the returned fraction R — 0.25. As 
done in Spitoni et al. (2010), those values were obtained by 
adopting the Salpeter (1955) IMF, and the stellar yields of 
Woosley & Weaver (1995) for oxygen at solar metallicity. 


3.1 Our model for a system formed by two 
isolated galaxies 

We consider the chemical evolution of the galaxyl with ini¬ 
tial mass Mtoti(O) = Mgi(O). We study here the effects of 
the interaction with another galaxy on the abundance of 
oxygen Zo,i(t) = Mo,i{t)/Mgi{t). 

We assume an environment dependent infall, originated 
by feedback episodes and galaxy-galaxy interactions. We as¬ 
sume that the infall is proportional to the gas outflows from 
nearby galaxies with wind parameter A. First we consider 
the case with the “enriched infall” proportional to the SFR 
of the galaxy2: 

Ax,lit) = eWx, 2 {t) = £A(1 - R)Mt), ( 9 ) 

where V' 2 (t) is the SFR of the galaxy2, and e is the frac¬ 
tion of outflowing gas from galaxy2 which reaches galaxyl. 
We assume in this work that e = 0.5. In fact, the outflows 
are generally bipolar and for a lobe approaching galaxy 1, 
the second lobe moves away from it. The value e=0.5 is an 
upper limit for the fraction of outflowing gas from galaxy2 
to galaxyl. Only if the outflow is well collimated along a 
narrow solid angle and only if galaxyl happens to intercept 
this galactic wind, e can be equal to 0.5. For example, from 
GALEX UV images of the starbust galaxy M81 given by 
Hoopes et al. (2005), it can be seen that the opening angle 
of the northern side wind is around 55 degrees. 

The Zo, 2 {t) metallicity of the infall is time and depends 
on the outflows of the galaxy2. 

The system of equations we have to solve for our galaxyl 
is the following : 

= (£A)(1 - R)Mt) 

= (1 - R)[e\'ip2{t) - ipi{t)) 

= (1 - R){Mt)[yo - Zo,iit)] +Mt)eXZo, 2 it)), 
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where is the SFR of the galaxyl. 

As mentioned above, Zo, 2 {t) is the metallicity of the 
galaxy2 which is suffering only gas outflows. We therefore 
need to solve another system of equations for the evolution 
of Zo, 2 {t) in galaxy 2 . 

Recalling Matteucci (2001) the system to be solved for 
galaxy 2 with only outflow of gas is: 


dMtnt2 

dt 

dMg2 

-- 

dMQ2 

dt 


= (1 — R)tp 2 {t) [ — Zo, 2 {t) + yo 


( 11 ) 

XZoAt)]- 


We recall that the outflow rate for the galaxy2 is W\^ 2 (t) = 
With equations (Hoj and (El) we study an isolated 
system formed by two galaxies when e = 1. In fact, in this 
case we have that = 0 . 

On the other hand, imposing e 7 ^ 1, we assume that part 
of the chemical enriched gas escaping from galaxy 2 , ends up 
in the in intergalactic medium (IGM). 


3.2 Model including primordial infall of gas and 
the interactions of 2 galaxies 

The majority of detailed chemical evolution models of galac¬ 
tic systems assumes that the galaxies formed by accretion of 
primordial gas from the IGM. Therefore, we examine here 
also the more realistic case where both galaxyl and galaxy 2 
suffer an inflow of primordial gas. We assume that the in¬ 
fall rate is proportional of the SFR of the galactic system 
as discussed in Section 2.1. Hence, for the galaxyl the pri¬ 
mordial infall rate is AA,i{t) = A(1 — R)'ipi{t), whereas is 
AA, 2 {t) = A(1 — R)'>jj 2 {t) for the galaxy2. The systems of 
equation we need to solve are the following ones: 


• galaxyl: 

' = (1 - R)ieXMt) + AV>i(f)) 

—^ = (1 - R){eX%l)2{t) -b (A - 

- R){i’i{t)[yo - Zo,i{t)] +'ip 2 {t)eXZo, 2 {t))- 

• galaxy 2 : 


dMfnt2 

dt 

dMg2 

dt 

dMo,2 

dt 


= {-X + A){l-R)ij2{t) 

(A-A-l)(l-R)V> 2 (t) 

= (1 - R)'ip 2 {t) [ — Zo, 2 {t) + yo 


(13) 

XZo,2{t)]- 


4 THE NEW ANALYTICAL SOLUTIONS 

In this Section we present our new analytical solutions. First, 
we discuss the case where the galaxyl suffers an enriched 
infall from an evolving companion galaxy (galaxy2). In the 
following we generalize this solution in the case of different 
SFEs for the galaxyl and galaxy2. We also present the case 
where a primordial infall is taken into account. The last re¬ 
sult is related to “environment” effects coupled with galactic 
fountains. 



Figure 1. Evolution of the metallicity Zq in terms of the gas 
mass fraction fj, for the closed-box model (blue solid line) and for 
a model with only gas outflow with A = 2 (magenta dotted line). 
Different horizontal solid green lines refer to different galactic 
time (1, 2, 4, 8, and 13 Gyr). 


4.1 Interactions with a nearby galaxy 

From the system of equations m we recover Zo, 2 {t), and 
following Matteucci (2001) we obtain: 

Zo,2it) = j^^rt[il + X)ft2"-X]. (14) 

Assuming a linear Schmidt (1959) law {ip = Sx Mg) (Pipino 
et al. 2014, Recchi et al, 2008), in the systems (IIOII and (Hill , 
we have the following expressions of the time evolution of 
the gas for the galaxyl (Mgi(t)) and the galaxy2 {Mg 2 {t))-. 
( 12 ) 

Mgi{t) = e-^^-^^^\Mgi{0)+eMg2{0) ), (15) 

Mg 2 {t) = Mg 2 ( 0 )e"A + l)(l-«)St^ 

Recalling that the gas fraction is defined as /r = 
Mg{t)/Mtot{t), we can write for the galaxy2: 

-1 _ Mtot2 Mtot2 

^ 1^ “ Mg2(0)e-(^+i)(i-«)^*' ^ ^ 

From the system El with the initial condition that 
Mtot 2 { 0 ) = Mg 2 ( 0 ), we have that: 

Mtot2{t) = ^Mg2{t) + Mg2{0) (1 - . ( 18 ) 

Finally, inserting eq. m in eq. El , we recover the following 
time dependent expression for /rj^(t): 

+ (M) 

Hence, the expression of Zo, 2 {t) as a function of the galactic 
time becomes: 

ZoAt) = In - A] = yoil-R)St.{20) 
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Figure 2. Left panel: The time evolution of the metallicity (oxy¬ 
gen) for the galaxyl when galaxy-galaxy interactions are taken 
into account and assuming an initial metallicities fixed at the 
value of Zo^i{0)=0.05 and Zo,2(0)=0. For all the cases we as¬ 
sume, A=0.4, e=0.5, S=1 Gyr“^. The magenta short dashed line 
represents the case with Mgi{0)/Mg2{0)=2, the green dotted line 
for Mpi(0)/Mg2(0)=l, the solid blue line represents the case with 
Mgi(0)/Mg2{0)=0.b. The long dashed grey line is the closed-box 
solution with the same initial metallicity. Right panel:-. Same of 
left panel but with an initial metallicity fixed at Zo,i(0)=0.1. 



Figure 3. The effect of the fraction of the outflowing gas (e) 
from galaxy2 to galaxyl on the time evolution of the metallicity 
(oxygen) for the galaxyl. We assume Zq^i{0)=0.1, Zq 2(0)=0, 
A=0.4, S=1 Gyr“^, and Mgi{0)/Mg2{0)=l. The blue solid line 
represents the model with e=0.5, the case with e=0.4 is reported 
with the dotted green line, the short dashed magenta line is the 
case with e=0.3. The models with e=0.2 and e=0.1 are in brown 
dashed dotted line and violet long dashed dotted line, respec¬ 
tively. The long dashed grey line is the closed-box solution with 
the same initial metallicity Zo^i(0)=0.1. 


It is worth noting that the evolution in time of Zo, 2 {t) 
does not depend on A in the specific case of the Schmidt 
(1959) SFR with /c=l, and is the same as in the closed-box 
one. In fact, for a closed-box model with the linear Schmidt 
(1959) law we have this expression for the gas mass: Mg(t) = 
. Following eq. Q we derive the expression 
of the metallicity Zcb as a function of time for the closed-box 
(c5). Therefore, we have that: 

Zcb{t) = yo In = j/o(l - R)St. (21) 

In Fig. m we draw the evolution of the metallicity Zq 
as function of the gas fraction y for both closed-box model 
and the one with linear outflows in the case of the linear 
Schmidt (1959) law, ip = SMg, with wind parameter A=2 
and S=1 Gyr“^. Horizontal lines indicate different galactic 
times during the galactic histories. 

Even if the closed-box model and the model with only 
outflow have the same metallicity in time, that same metal¬ 
licity is reached for different values of the gas fraction y. 

The metallicity Zo, 2 {t) expressed in eq. (1201) can 
be inserted in the system la, and recalling that 
~ , the third equation of the system (HOI) 

can be rewritten as: 

= (1 - R)x 

(ipi{t)yo + eip2{t)[- XZo,i{t) + XZo,2{t)])). ( 22 ) 


Hence, we assume that the SFR follows a Schmidt (1959) 
law: ipi = S X Mgi and ip 2 = S x Mg2, and the differential 
equation we need to solve when we include eq. ([20)) in eq. 
(021 is: 


dZo,i{t) 

dt 


= {l-R)Sx 


{vo + [ - eXZoAt) + eAj/o(l - R)St]y 


(23) 


The final expression for Zo,i{t) in a case where the 
galaxyl is affected by the interaction with another galactic 
system (galaxy2), and with Zq,2(0) = 0 is the following one: 


Zo,i(t) = Syo{l — R)t+ ■ 


Z 0 , 1 ( 0 ) 


1 + e 


Mg2(0) 
Mgl(O) ' 


I _ g-A(l-fl)St 


-.(24) 


We see that the closed-box solution is still recovered 
when the initial metallicity Zq, i(0) for the galaxyl is con¬ 
sidered equal to zero. The reason for this behaviour is that 
the outflow from galaxy2 has the same metallicity of the 
whole galaxy, i.e. it is expressed by eq. (20). Since the two 
expressions (eqs. 20 and 21) are identical, and since metallic¬ 
ity is not an additive quantitiy, the expression Zo,i{t) in eq. 
(24) is expected. In Fig. Owe show the time evolution of the 
metallicity (oxygen) where we assume an initial metallicity 
different from zero in the galaxyl and as function of differ¬ 
ent Mgi{0)/Mg2{0) ratios: 2, 1 and 0.5. We considering two 
different initial metallicities Zo.i(0)=0.05 (left panel) and 
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Zo,i{0)=0.1 (right panel). We are aware that these values 
are extremely large, but here we only like to show the trends 
of considering different pre-enriched values. We see that at 
early times the dilution effect is prominent because the in¬ 
falling gas has a much lower metallicity of that galaxyl, and 
as expected the lower the ratio Mgi( 0 )/Mg 2 ( 0 ) is, more the 
dilution effect is important. 

In Fig. [ 3 ] the effects of the fraction of outflowing gas 
e on the time evolution of the galaxyl oxygen abundance 
obtained by eq. ([24]) are shown. Assuming Zo,i(t)=0.1, 
Mgi(0)/Mg2(0) = l, A=0.4, and S= 1 Gyr“^, we present 
different models varying e: e = 0.5, 0.4, 0.3, 0.2, 0.1. As ex¬ 
pected, the dilution effect is larger assuming higher e values. 

We discuss now the case where the initial metallicity 
for galaxyl and galaxy2 are Zo,i{0)=0, Zo, 2 ( 0 ) 7 ^ 0, re¬ 
spectively. In this case the differential equation we need to 
solve is: 

dt ^ ’ 


,yo+ 


Mg2(t) 


Mgl(t) L 


-eA.Zo,i(t)-t-eA^.Zo,2(0)-t-yo(l— R)St^ ^ -(25) 


The new solution is: 


Zo,i(t) = Syo{i—R)t+ 


Mgl(0) 


1 _ g-A(l-K)S£ 


1 + ^^S2(.0) A _ -A(l-H)StA 


.(26) 


In Fig. |4] we show the evolution of the oxygen abundance 
of galaxyl using the new analytical solution presented in 
eq. (I2nj with A=0.4, £=0.5, and Zo,i(0)=0, considering two 
different initial metallicities Zo,2(0)=0.05 and Zo,2(0)=0.1. 
In this case we use such high values for the pre-enrichment 
in order to better visualize the general trends and the effects 
of different Mgi(0)/Mg2(0) ratios. As expected, at variance 
with Fig. [2] the smaller is the Mgi (0)/Mg2(0) ratios, the less 
efficient the chemical evolution is for galaxyl. 

Overall, the general solution for the evolution of the 
oxygen abundance for the galaxyl with both initial metal¬ 
licities Zo, 1 ( 0 ), Zo, 2 ( 0 ) different from zero, is given by the 
following expression: 


Zo,i(t) = Syo{i — R)t + (27) 

Zo,i{0) + Zo.2{0)e^^ (1 - 

The general solution reported in eq. (|27ll describes the more 
realistic scenario where galaxies exhibit different chemical 
enrichment rates and therefore they posses different metal¬ 
licity when the interaction starts. It seems plausible that 
galaxy2 presents a faster chemical enrichment in systems as 
the M81 group, where large and metal rich galaxies (M81, 
M82, and NGC3077) interact and eject gas. In other envi¬ 
ronments, it can be more reasonable that only dwarf galaxies 
would show galactic winds, and therefore in this situation 
galaxy2 is less chemically enriched than galaxyl. In con¬ 
clusion, both cases Zo,2{0)>Zo,i{0), and Zo, 2 ( 0 )<Zo,i{0) 
have physical meaning. Recently, in Recchi et al. (2015) it 
was shown that the values of initial metallicities depend on 
the level of the interaction between galaxies. 




Figure 4. Left panel: The time evolution of the metallicity (oxy¬ 
gen) for the galaxyl when galaxy-galaxy interactions are taken 
into account and assuming an initial metallicities fixed at the 
values of Zo,i{0)=0 and Zo,2(0)=0.05 . For all the cases we 
assume, A=0.4, e=0.5, S=1 Gyr“^. The magenta short dashed 
line represents the case with Mgi(0)/Mg2(0)=2, the green dot¬ 
ted line for A4gi(0)/Mg2(0)=l, the solid blue line represents the 
case with A4gi(0)/Mg2(0)=0.5. The long dashed grey line is the 
closed-box solution with the same initial metallicity. Right panel:-. 
Same of left panel but with the initial metallicity fo galaxy2 fixed 
at ^o,2(t)=0.1. 


4.2 Some applications to real cases 

The evolution of the stellar mass content in the galaxyl 
M*p(t), can be simply inferred by the relation M*,i(t) = 
Mtoti{i) — Mgi(t). The time evolution of the total mass of 
galaxyl is given by the first equation of system 

(HQi, and with the initial condition Mtoti{ 0 ) = Mgi(O) we 
have that: 

Mtotiit) = Mgi(O) + _ ^-Hi-R)sty (38) 

In Fig. [ 5 ] we compare the time evolution of 
log{0/H)+12 as function of stellar mass for the closed-box 
model with a model where environment effects are taken 
into account. The \og{0/H) + 12 quantity is recovered from 

Zo,i{t) using the following expression: \og(^Zo,i{t)/{lQ x 

0.75)^ +12, where 0.75 is the assumed fraction of hydrogen. 

The assumed initial mass in each model is 10^'^ M©. In the 
last case, at a fixed stellar mass value the system shows a 
smaller oxygen abundance compared to the closed-box evo¬ 
lution. 

In Fig. [6] we study the effects of the environment on 
the mass-metallicity (MZ) relation by comparing with the 
observed relation of Kewley & Ellison (2008) for Sloan Digit 
Sky Survey (SDSS) star-forming galaxies. To estimate the 
amount of gas that resides in each star-forming galaxy, and 
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Figure 5. The evolution in the framework of the closed-box 
model of the abundance \og{0/H)-\-12 as a function of the stellar 
mass M* for the closed-box model for a galaxy with total mass 
equal to 10^® Mq is drawn with the dotted green line. With the 
solid blue line we show the results for the galaxyl with a initial 
mass of 10^^ when we take into account the enriched 

infall from a companion galaxy with the same initial mass and 
the wind parameter equal to A=2, and e=0.5. The initial oxygen 
abundances are ^o,i(0)=-^O,2(0)=0. 

the gas fraction ^ as a function of the galactic stellar mass, 
we use the method described in Spitoni et al. (2010). 

We determine the cold gas mass of each galaxy on the 
basis of its SFR, using the following inverted Kennicutt 
(1998) relation, which links the gas surface density to the 
SFR per unit area: 

- (ssliiFr)"'™ 

where the gas density is expressed in M 0 pc“^, and 
the SFR surface density E* in M© yr“^ kpc~ . The gas mass 
Mgas (in Mq) is given by: 

— Egtas X 27vRI (30) 

where Rd is the scaling radius calculated as in Mo et al. 
(1998). At this point we have a relation between ^ and the 
stellar mass for each considered galaxy. 

First, we follow the evolution of the closed box mod¬ 
els with zero initial metallicity for different initial masses 
Mtoti0)=Mg{0)=Mtotit)=8.8xl0\ 5x10®, lO^®, 5xl0“, 
and 10^^ M 0 , respectively. For each model we compute the 
time evolution of Zo{t), /r(t), and Thus, we con¬ 

sider the time Imz where the stellar mass Mt,{tMz) and 
the corresponding jj, belong to the fit reported in Fig. 2 
of Spitoni et al. (2010). At this point the metallicity at 
the time tMZ is computed using eq. (21). In Fig. [ 6 ] we 
show the closed-box model results for systems with initial 
masses Mtot{0) = Mg{0) = = 8.8x10®, 5x10®, 10“, 

5x10^®, and 10^^ Mq, respectively. 


Figure 6. The observed mass-metallicity relation and related 
standard deviation of Kewley & Ellison (2008) are indicated with 
the dotted grey lines. With the green circles we show the closed- 
box results for galaxies with total masses: 8.8x10®, 5x10®, 10^®, 
5x10^®, 10*^^ A/q, respectively. We compare them with model re¬ 
sults where we consider the interactions with a companion galaxy 
with the same initial masses. For all the models it is assumed that 
the initial oxygen abundances are Zo,i{0)=Zo,2(0)=0. With the 
magenta hexagon models with the wind parameter fixed at the 
value of A = 0.1 are presented. The cases with A = 0.5, A = 2, 
and A = 6 are represented with blue pentagons, red squares, and 
brown triangles respectively. 


The closed box model is not able to reproduce the ob¬ 
served MZ relation, and in Spitoni et al. (2010) we concluded 
that a galactic wind rate increasing with decreasing galac¬ 
tic mass or a variable IMF are both viable solutions for 
reproducing the MZ relation. We consider also galactic sys¬ 
tems with the same initial masses of the closed-box models, 
but taking into account environment effects: an enriched in¬ 
fall from a galaxy with the same mass and with the same 
SFE with zero initial metallicities for galaxyl and galaxy2: 
Zo,i{0)=Zo, 2 ( 0 )—0. 

In Fig. [ 6 ] we show the results when we consider four 
different inflow parameters: A=0.1, 0.5, 2, 6 . We see that at a 
fixed stellar mass the metallicity of the galactic system drops 
down in presence of inflow. We are in agreement with the 
work of Torrey et al. (2012), where the effects of dynamical 
interaction of companion galaxies were studied. In Spitoni 
et al. ( 2010 ) we have already proved that a combination of 
primordial infall and outflows with variable wind parameters 
as functions of the stellar mass is required to reproduce it. 
Here, we show the effects of this interaction compared to the 
closed-box case. 

In Torrey et al. (2012) galaxy-galaxy interactions do not 
change the slope of the MZ, but the MZ is just shifted down 
in parallel at smaller metallicities. In our case, if we con¬ 
sider a constant A (in Fig. [ 6 ] connecting the same type and 
colour points) for all the galaxy masses we have a fiatten- 


© xxxx RAS, MNRAS 000, Q]-?? 







8 E. Spitoni 


9.2 


—I-1-1-1-1-1-1-1 1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1-1— 


9 - 


w 8.8 

I 

+ 

ac 
\ 
o 


8.6 - 


8.4 


8.2 - 


t 


f 

* 




I I I I I I I I I I I I I I I I I I I I I I I I I I 


f 


closed box 

X-0.1 

X-0.6 

X-6 

_L 


8.5 9 


3.5 10 

log(M./Mg) 


10.5 11 



Figure 7. The observed mass-metallicity relation and related 
standard deviation of Kewley Sz Ellison (2008) are indicated with 
the dotted grey lines. With the green circles we show the closed- 
box with the variable yo values of Spitoni et al. (2010) for galaxies 
with total masses as in Fig. [6] We compare them with model re¬ 
sults where we consider the interactions with a companion galaxy 
with the same initial masses. For all the models it is assumed that 
the initial oxygen abundances are Zo,i(^)=Zo,2(^)=0. Models 
with different wind parameters A are indicated with symbols as 

Fig. [el 
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Figure 8. The observed mass-metallicity relation and related 
standard deviation of Kewley &: Ellison (2008) are indicated with 
the dotted grey lines. Models are indicated with symbols as Fig. 
[e] assuming that the initial oxygen abundances are for galaxy 1 
^O,i(0)=05 for galaxy2 Zq^ 2(0)= 5 x 10~^, respectively. 
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Figure 9. We show the metallicity Zq ^ in the case of enriched 
infall from the evolution of galaxy2 as a function of ln(/i“^). 
We fix 6=0.5, A=2, and consider different ratio Mgi(0)/Mg2(0) 
values: 1 (green dotted line), 0.1 (magenta short dashed line), 
0.01 (grey long dashed line). We compare our results with closed 
model case (solid blue line). The initial oxygen abundances for 
the galaxyl and galaxy2 are Zq i(0)=Zo 2(0)=0. 


ing of the MZ slope compared to the closed-box. However, 
from Spitoni et al. (2010) we know that more massive galac¬ 
tic systems are associated to smaller outflow episodes and 
therefore smaller wind parameters. Hence, if we consider a 
variable A (decreasing towards higher mass galaxies) we pre¬ 
serve the slope of the MZ relation. 

Our results are also in agreement with the observed 
depression of 0.05 -0.10 dex found by Ellison et al. (2008) 
for the MZ relation in interacting galaxies, if we consider a 
variable wind parameter A with a maximum value smaller 
than A=6. 

In Fig. [7] we show the MZ relation obtained by the 
closed-box models with the same prescriptions adopted in 
Fig. [6] but with the variable yields ycLj of Spitoni et al. 
(2010) which is able to perfectly reproduce the MZ relation. 
Again, the fact of considering “environment” effects leads 
to a drop in the metallicity of the galaxyl at a fixed stellar 
mass compared to the closed-box models. The fit with the 
observations deviates for A > 0.1, and this might be a hint 
that the kind of interactions treated in this paper cannot be 
very general, but it surely applies for specific systems (such 
as the M81 group). 

Finally, in Fig. [S] we show the effects of different ini¬ 
tial metallicities Zo, 2 {Q) of galaxy2 on the MZ relation. 
Assuming the same model parameters of Fig. but with 
•^0,2(0) = 5 X 10“®, we see that at a fixed stellar mass, the 


^ We remind that varying the yield for stellar generation corre¬ 
sponds to vary either the stellar nucleosynthesis or the IMF. 
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decrease of the metallicity is less prominent compared Fig. 
[6] especially for small stellar masses. 

Edmunds (1990) studied the effects of gas flows on the 
chemical evolution of galaxies showing that it does exist a 
forbidden area in the plane log(/r“^) versus metallicity above 
the simple model solution, for models with inflow of unen¬ 
riched gas or outflows. Here, we want to test if this result 
is still valid in the case of a time dependent enriched in¬ 
fall from a companion galaxy with initial metallicities for 
galaxyl and galaxy2 set to zero: Zo,2{0) = Zo,i{0)=0. The 
expression of /.ii(t) as a function of time is: 


Mi(f) = 




1-b 


Ae M„2(0) 
1-t-A Mgi(O) 


(^1 _ 


(31) 


We need an expression for /ii(t) as a function of the metal¬ 
licity Zo- We have seen that the solution Zo,i{t) in the case 
of zero initial metallicity is (1 — R)Syot. Therefore, we can 
rewrite eq. la as: 


ni{Zi) 




1-b 


Ae M„2(0) 
1-l-A Mgi(O) 




(32) 


In Fig. [9] we show the metallicity Zo,i in case of en¬ 
riched infall from the evolution of galaxy2 as the function 
of ln(^“^). We fix e=0.5, A=2, and consider different ra¬ 
tio Mgi{0)/Mg2{0) values:!, 0.1, 0.01. We see that all the 
studied cases are in the forbidden area defined by Edmunds 
(1990), and they evolve similarly to the closed box model 
with small differences. This result, at variance with Ed¬ 
munds (1990) is due to the fact that here is the first time in 
which these particular cases are studied and the theorems 
T(2) and T(3) of Edmunds (1990) cannot be applied because 
of the presence of enriched infall. 


4.3 The analytical solution including both 

galactic fountain and environment effects 

In Recchi et al. (2008), new analytical solutions in the frame¬ 
work of differential winds were presented, namely galactic 
winds in which the metals are ejected out of the parent 
galaxy more efficiently than the other elements. The exis¬ 
tence of differential winds has been first introduced in the 
context of chemical evolution of galaxies by Pilyugin (1993) 
and Marconi et al. (1994). A realistic case of variable in¬ 
fall metallicity is represented by the situation in which the 
metallicity of the infalling gas is set to be always equal to 
the one of the galactic wind. This condition implies that the 
very same gas that has been driven out of the galaxy by en¬ 
ergetic events can subsequently fall back to the galaxy, due 
to the gravitational potential well. 

This kind of duty cycle is called galactic fountain 
(Shapiro & Field 1976, Bregman 1980). Spitoni et al. (2008, 
2009) showed the effect of galactic fountains on a detailed 
chemical evolution model (where the instantaneous recycling 
and mixing approximation were relaxed) of the Milky Way. 
They discussed the delay in the chemical enrichment due to 
the fact that the gas takes a finite time to orbit around the 
Galaxy an fall back into the disk. In Recchi et al. (2008) 
an analytical solution was presented in the case for galactic 
fountains. 

Here, we want to generalize the results of Recchi et al. 



Figure 10. In this sketch the gas flows related the galaxyl are 
represented when galactic fountain are considered. The enriched 
infall with the parameter A is drawn with the yellow arrow. The 
outflowing gas from galaxyl with parameter lo and the fraction 
which fall back (SI) are in gray arrows. Our model is in the frame¬ 
work of differential wind theory when metals are more easily chan¬ 
neled out. 



Figure 11. We show the time evolution of the metallicity Zo,i in 
presence of galactic fountains when mass ratios Mgi{0)/Mg2{0) = 
1 and SFE 5=1 Gyr“^, using the new analytical solution pre¬ 
sented in eq. IIAlI l with e=0.5, and A=0.4, S! = 0.1 and oj = 0.2. 
The initial oxygen abundances for the galaxyl and galaxy2 are 
Zo,l{0)=Zo,2{0)=0- The model with a= 4.5 is drawn with the 
blue long solid line, the model with a= 3 is represented by dotted 
green line solid line, the case with a= 2 is drawn with with the 
short dashed magenta line. With the cyan long dashed line the 
model with a= 1.5 is represented. 


(2008), and find a new analytical solution for the metallicity 
(oxygen) in the galaxyl when the galactic fountain effect 
is included together with the the enriched infall from the 
galaxy2. To consider galactic fountains we have to take into 
account an outflow episode and a new infall one caused by 
SN explosion events in the galaxy, with parameters ui and 
Q. The metallicities of the outflows an fountains are: 


Zui = Zo,iQ:a;(l — RNi, 

(33) 

Zn = Zo,io:itl{l — R)'ipi, 

(34) 
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Figure 12. We show the time evolution of the metallicity Zq^i 
when we consider different SFEs S 2 for the galaxy2 following the 
new analytical solution presented in eq. goj. We consider the 
case with Mgi(0)/Mg2(0)=l, e=0.5, and A=0.1. The SFE for the 
galaxyl is fixed at the value of 51 = 1 Gyr“^ . The model with 
52 = 1 Gyr~^ is drawn with the blue solid line, the model with 
52=2 Gyr“^ is represented by the magenta short dashed line, 
With the long dashed grey line the model with 52=0.5 Gyr“^ 
is represented. Finally, the model with 52=0.1 Gyr“^ is drawn 
with the dotted green line. The initial oxygen abundances for the 
galaxyl and galaxy2 are Zq i{0)=Zq 2(0)=0. 



Figure 13. We show the time evolution of the metallicity 
Zo,i when we consider different mass ratios Mgi(0)/Mg2{0), 
whereas SFEs 5i and S 2 have been fixed at the values of 1 
and 0.1 Gyr”^, respectively. We use the new analytical solution 
presented in eq. Hot with e=0.5, and A=0.1. The model with 
'^9l(0)/'^92(0)= 10® is drawn with the grey long solid line, the 
model with Mgi{ 0 )/Mg 2 ( 0 )= 10® is represented by the blue solid 
line, the case with Mgi(0)/Mg2{0)= 1 is drawn with with the 
dotted green line. With the magenta short dashed line the model 
with Mgi{0)/Mg2{0)= 10“^ is represented. The initial oxygen 
abundances for the galaxyl and galaxy2 are Zq i{0)=Zo 2(0)=0. 


We are studying galactic fountain effects in the framework 
of differential winds as done in Recchi et al. (2008) where it 
was introduced the parameter a > 1 , in order to take into 
account for the fact that metals are more easily channeled 
out from the parent galaxy compared to the unprocessed 
gas. 

In Fig. [To] we present a sketch of the gas flow patterns 
for the galaxyl when galactic fountain are considered. It is 
shown the gas inflow from the galaxy 2 with the associated 
wind parameter A, the outflow from galaxyl (tj), and the 
fraction which comes back (Q). 

Melioli et al. (2015) using three-dimensional hydrody- 
namical simulations investigated the impact of SN feedback 
in gas rich dwarf galaxies and the formation of galactic foun¬ 
tain and outflows. They found a similar circulation flow of 
the one presented in Fig. 1101 galactic fountains is generally 
established and the metal-rich SN ejecta is instead more sig¬ 
nificant fraction (25-80 percent) is vented in the intergalactic 
medium, even. 

It must be underlined that the treatment of the galactic 
fountains is actually analogous to what is done in Recchi 
et al. (2008) with oj replacing A and fl replacing A. The 
difference is in the role of galaxy2. 

The system of equations we have to solve for the galaxyl 
is the following: 


[ = (1 - + {^- 
j = (1 - + {n-uj- i)V'i(i)) 

[ ^{i-R){Mt)[yo + QZoAt)] +Mt)e\Zo,2{t)). 

with 0 = —1 -b aQ. — auj. 

The differential equation of the metallicity in terms of 
oxygen for the system is: 


dZo,i{t) 

dt 


{l-R)Sx 


(j/o+^o.i(i)x+^^[-eAZo,i(t)+cAyo(l-R)5t]),(36) 


with X = (a — l)(fl — u;). It can be easily proved that the 
evolution in time of the mass of gas of the galaxyl is given 
by: 


Mglit) 


1 “h ^ 


(^A(l-e“^^+^"“'^^“^>®‘)eMg2(0)-bM9i(0)(A-bfl-a;)).(37) 

In eq. CD) of the Appendix A we show the complete ex¬ 
pression of the analytical solution in presence of galactic 
fountains. 

In Fig. [11] the effects of different values of the param¬ 
eter a are tested. The higher a is, the less the galaxy gets 
enriched. This is because we assumed that fl < cu and con- 
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sequently the metals (in this case oxygen) which escape, are 
larger than the ones which rain back into the galaxy. 


4.4 Model with different star formation 
efficiencies 


Finally, we test the case where the interacting galaxies have 
different SFEs. It is well known that different galactic sys¬ 
tems can be characterized by different SFEs (Matteucci 
2001), and generally, higher SFEs can be associated to more 
massive systems. Here, we show the new analytical solution 
for the galaxyl with a SFE receiving enriched outflow of 
gas from another galaxy associated with a SFE S 2 - There¬ 
fore, with the new SFRs tpi = Si x Mgi and 'ip 2 = S2 x Mg2 
the eq. m becomes: 


Mgl{t) 


g-(l-H)Sii 

5i - (A -h 1)S'2 ^ 


|Mgi(0)(Si-(A-hl)S2) -f eAMg2(0)S2 (e^* - l)|, (38) 


with 5 = (1 — R){Si — (A -I- 1 ) 52 ). It is trivial to see that 
when Si = S 2 eq. (I38II is identical to eq. (I15II . The differential 
equation we have to solve to have the metallicity evolution 
of the galaxyl in presence of different SFEs is: 


dZo,i{t) 

dt 


{1-R)x 



+ [ - eXZo,lit) + eAyo(l - R)S 2 f]) . (39) 


Finally, the new analytical solution in this case is the fol¬ 
lowing: 


Zo,iit) 


1 

(Si - (1 -h A)^) ^ 



e\S2yo 

A/f) 

Mg2(0) 

B{t) 


( l-fe )e-^>S'2-b 

'si - (1 + A)5'2) 


(40) 


where A{t) and B{t) are the following time dependent terms: 


A{t) = S2 + e^‘|si - S 2 -h (-1 + i?)S2( - Si -h S2(l -f A))f| + 
-h Si|-1 -f (-1 -f i?)(Si -{1 + A)S2)f|, (41) 

B{t) = (-Si -h S 2 -f AS 2 )" ((-1 -f R)Sityo - Zo,liO)) ■ (42) 

We note that the result depends on the ratio between 
the initial masses of galaxyl and galaxy2. Moreover, as ex¬ 
pected, when we impose Si = S 2 in eq. (|40ll. the resulting 
Zo,i{t) is identical to the one obtained with eq. (I24II . In Fig. 
na we test the effect of different SFEs for the galaxyl and 
galaxy2 when we consider the same initial gas mass ratio 
(Mgi(0)/Mg2(0)= 1). We tested different values for S 2 : 2, 1, 

0.5, 0.1 Gyr“^. For the galaxyl we assume Si = l Gyr“^ in 
all models. 

We obtain, as expected, that smaller values of S 2 lead to 
a smaller chemical enrichment of the system. In the case of 
S 2 = 0.1 Gyr“^ the metallicity of the system even decreases 
at later galactic times due to the strong dilution effect of 



Figure 14. We report the time evolution of the metallicity Zq 1 
when we take into account also a primordial infall with wind pa¬ 
rameter A. We consider the case when Mgi(0)/Mg2(0)=l, e=0.5, 
A=0.3, and S=1 Gyr~^. The initial oxygen abundances for the 
galaxyl and galaxy2 are Zq i{0)=Zq 2(0)=0. The model with 
A=0.5 is drawn with the long dashed grey line, the model with 
A=0.7 is represented by the magenta short dashed line, the case 
with A=0.9 is drawn with the dotted green line. With the solid 
blue line the model with A=1.5 is represented. 


the infalling gas. In fact, we can consider it like pristine 
gas because of, as expected, the small SFE value of galaxy2 
compared to galaxyl. 

In Fig. I13[ fixing and S 2 at the constant values of 
1 and 0.1 Gyr“^, respectively, we show the effect of differ¬ 
ent initial gas mass ratios for the galaxyl and galaxy2. We 
consider four cases: Mgi( 0 )/Mg 2 ( 0 )= 10^, 10^, 1, and 10“^. 
Because of the choice of an extremely low SFE for galaxy2 we 
expect that the most important deviation from the closed- 
box solution is obtained with larger Mg 2 { 0 ) values, hence 
smaller Mgi(0)/Mg2(0) ratios. In Fig. I13l we confirm it, and 
the model with Mgi{ 0 )/Mg 2 { 0 )— 10^, is almost identical to 
the closed-box evolution model. 


4.5 The effects of the primordial infall 


We present here the results in the presence of a primordial 
infall of gas for both galaxyl and galaxy2. First, we need to 
compute the time evolution for the metallicity Zo, 2 (t) when 
a primordial infall of gas is considered. To this aim we need 
the expression of the mass fraction y. 2 {t) as a function of 
time. From the system GSl), with the condition that the total 
initial mass is equal to the gas content Mtot 2 { 0 ) = Mg 2 { 0 ) 
we have that: 

^‘°«(t) = j^^j^Mg2{t) - Mg2(0) ^_^_^ . (43) 


The mass fraction p 2 = Mg 2 {t)/Mtot 2 {t) can be written as: 


U 2 




A-A 




A-A-1 A-A-1 


) ■ (44) 
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Figure 15. We show the time evolution of the metallicity Zq^i 
when we consider different SFEs S 2 for the galaxy2 coupled with 
a primordial infall of gas following the new analytical solution 
presented in eq. llA^ of Appendix A. We consider the case when 
Afgi(0)/Mg2(0)=l, and e=0.5, and the SFE fro the galaxyl is 
fixed at the value of 51 = 1 Gyr~^ . All the models assume A=0.3 
and A=0.7. The initial oxygen abundances for the galaxyl and 
galaxy2 are Zq^i(0)=Zo,2(0)=0. The model with 52=2 Gyr~^ 
is represented by the magenta dotted line, the model with 52=1 
Gyr“^ is drawn with the blue solid line, the model with 52=0.2 
Gyr“^ is drawn with the green short dashed line. Finally, the 
model with 52=0.1 Gyr”^ is drawn with the grey long dashed 
line. 

We recall that the general solution as function of the gas 
fraction for a system with primordial infall and outflow is 
given by eq. © with Za = 0. Finally, if we insert eq. (|44|1 
in eq. di), we obtain that 

Zo ,2 = ^ (1 - . (45) 

The new expressions for the time evolution of Mg\{t) and 
Mg 2 {t) are, respectively: 

Mgi(f) = |, (46) 

Mg2{t) = (47) 

The differential equation we should to solve in presence 
of a primordial infall, is the following one: 

= (1 - R)S{yo - RZoAt)) + 

-c^o.i(t)A|.(48) 

The solution of eq. (1481) for the time evolution of the oxygen 
abundance for the galaxyl is then given by: 

Zo,iA) = —7- 

K\Mgi(S)) + (1 - e^(-i+«)S'*)£Mg2(0) j 


[A2{t)+ [Mgi{Q) + eMg2{Q)]yo-B2{t) + C2{t)^, (49) 

with 

A2{t) = (-l + e^^"^+^^®‘)e^'"^+'^'®*eMg2(0)yo, 

B2{t) = e(^)'-i+«)®‘eMg2(0)y„, 

and 

C2(t) = [yo - A^o,i(0)]. 

It is easy to show that it is possible to recover the solution of 
eq. (1241) presented in Section 4.1 for limA^-o Zo,i{t). In Fig. 
na we show the time evolution of the metallicity Zo,i when 
Mgi{0)/Mg2{0) = l, e=0.5, and the SFE for the galaxyl and 
galaxy2 are fixed at the value S=1 Gyr“^ . Assuming for all 
the models A=0.3, we test different values for the the infall 
parameter A which is associated with the primordial infall: 
0.5, 0.7, 0.9, 1.5. 

We note that the effect of dilution of a primordial gas 
infall overwhelms the enriched infall from the companion 
galaxy, and the time evolution of the metallicity for the 
galaxyl deviates substantially from the closed-box solution 
(we have proved in previous Sections that the evolution in 
time of a system with only enriched gas from a compan¬ 
ion galaxy follows the closed-box solution) even with a pri¬ 
mordial infall parameter comparable to the one of enriched 
infall. 

We also consider the case with different SFEs for the 
galaxyl and galaxy2 in presence of primordial infall. We do 
not enter into details about the procedures and the relations 
used to recover the new analytical solution Zo,i{t) , and we 
only show the solution in ea. (IA2l) of Appendix A. 

In Fig. [TS] we show the effects of different values for 
SFEs of the galaxy2 on the chemical evolution of galaxyl in 
presence of primordial infall. We consider the case when the 
mass ratio is Mgi(0)/Mg2(0) = l, and e=0.5, and the SFE 
for the galaxyl is fixed at the value of = 1 Gyr~^ . All 
the models assume A=0.3 and A=0.7. We consider different 
models with different S 2 values: 2,1, 0.2, 0.1 Gyr“^. In Fig. 
I15l it is clearly shown that all the models with different SFEs 
show more or less the same time evolution for the Zo,i- 

Comparing the models with 52=0.2 Gyr“^ and 52=0.1 
Gyr“^ we note that at late times, as expected, the chemical 
evolution of galaxyl, which suffers the enriched infall from 
galaxy2 with 52=0.2, shows a higher metallicity. This be¬ 
havior is inverted at early times. This is due to the way in 
which we consider the primordial infall. In fact as shown in 
eq. (5), the primordial infall is assumed to be proportional to 
the SFR of the galaxy, therefore at early times the galaxy2 
with 52=0.2 Gyr”*^ suffers a larger dilution effect than the 
same system with 52=0.1 Gyr“^. 

We conclude that in presence of primordial infall the 
chemically enriched gas coming from galaxy2 has not a big 
effect on the chemical evolution of the galaxyl. 

In this work we consider only primordial infall from the 
IGM, whereas in Peng & Maiolino (2014b) it was studied 
how global environment properties (overdensity of galax¬ 
ies) can modify the metallicity of the IGM. The observed 
strong correlation between over-density and metallicity for 
star-forming satellites suggests that the gas infall is getting 
progressively more metal-enriched in dense regions. 
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5 CONCLUSIONS 

In this paper we presented new analytical solutions for the 
evolution of the metallicity (oxygen) of a galaxy in pres¬ 
ence of “environment” effects coupled with galactic foun¬ 
tains, and primordial infall of gas. The main results of our 
study are the following ones: 

• If we consider a linear Schmidt (1959) law for the star 
formation rate we have the same time evolution of the metal¬ 
licity for both the closed-box model and the one with only 
outflow. Therefore, in the last case the result does not de¬ 
pend on the wind parameter A. This result is holding when 
the outflow is not differential. It is worth to note that even if 
at the same time those two systems show the same content 
of metals they show different gas fractions as shown in Fig. 

m 

• The new analytical solution for the evolution of the 
galaxyl where we consider the enriched inflow of gas from 
the galaxy2 with initial metallicities Zo,i{0) (galaxyl) and 
^ 0 , 2 ( 0 ) (galaxy2) is: 

Zo,iiy) = 'S'i/o(l “ R)t+ 

ZoAO) + ^o, 2 ( 0 )e^^ (1 - 

where: A is the wind parameter, e is the fraction of the out¬ 
flowing gas from galaxy2 to galaxyl, S is the star formation 
efficiency, and is the ratio between the initial gas 

masses of the two galaxies. 

• When we consider the time evolution of a galaxy in¬ 
cluding an enriched infall due to the interactions with a 
nearby galaxy assuming no pre-enriched gas for both galax¬ 
ies, the chemical enrichemnt is less efficient than in the 
closed-box solution at a fixed stellar mass. Moreover, the 
mass-metallicity relation for galaxies which suffer a gas in¬ 
fall from an evolving galaxy, is shifted down at lower Z val¬ 
ues if compared to the closed-box model results for isolated 
galaxies. 

• We show a new analytical solution where we consider 
galactic fountain effects in the framework of differential 
winds coupled with the interaction with a nearby galaxy. 

• We presented the new solution in case of different SFEs. 
If the infall gas originated by a nearby galaxy has a smaller 
SFE than the galaxy suffering the infall, the recipient galaxy 
will show a less efficient chemical enrichment. In this case we 
showed that the smaller is the ratio between the initial gas 
masses of the recipient and donor galaxy, the less efficient 
the chemical evolution is. 

• A new analytical solution when a primordial infall is 
coupled with interactions with a nearby galaxy is presented. 
The effect of dilution of a primordial gas infall overwhelms 
the enriched infall from the chemical evolution of a compan¬ 
ion galaxy even with a primordial infall parameter compa¬ 
rable to the one of enriched infall. We have also shown how 
different SFEs for the companion galaxy does not affect the 
chemical evolution of galaxyl in presence of primordial in¬ 
fall. 
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APPENDIX A: NEW ANALYTICAL SOLUTIONS WITH GALACTIC FOUNTAINS AND 
PRIMORDIAL INFALL 


• We report here the analytical solution for the time evolution of the oxygen abundance Zo,i when galactic fountain and 
the enriched infall of gas from a companion galaxy are taken into account. This solution is obtained when the initial metallicity 
of the galaxyl is assumed equal to zero. 


(-l + a)(A + a(n-cu))2A3(t)(n-cu) ^ 

(a'^il - + X^B3{t) + A^(n - uj)C3{t) + A(fl - ujfD3{t)^ 

with 

A3{t) = A(Mgi(0) + (1 - + Mgi{Q){Q. - to) 

B3{t) = (l - (M9i(0) + eMg2(0)) + (-1 + a)e^^+“'“““”'"^+^^®*£Mg2(0)(n - a;)(-l + R)St 

C3{t) = (l - (^Mgi{0) + 2aMgi{0) + 2aeMg2{Q)^ + (-1 + - a;)(-l + R)St, 

D3{t) = -eMg2(0) + aE3{t) + (-1 + a)e^^+“^^"‘"''^“^+^^®‘eMg2(0) - 1 + a + a(fl - a;)(-l + R)St^, 

E3{t) = (2 + a)Mgi(O) + 2eMg2(0) - )(-i+fl)st^(2 + a)Mgi(O) + aeMg 2 ( 


• The following solution is related to the case with a primordial infall coupled with an enriched infall from a companion 
galaxy, considering different SFEs: 

A4(t) + B4(t) 


Zo,lit) = 


with: 


A(S'i - (1 + A)52)(1 - e'J‘)eAMg2(0)S'2 + Mgi(0)((-1 + A)Si + (1 - A + A)S'2) 


(A2) 


q={-l + R) ((-1 + A)^! + (1 - A + A)S' 2 ) , 

q2 = q/(-l + R), 

A* = A(-l + i?)S'i, 
e\Mg2iO)S2, 

A4it) = eAMg2(0)S'2([l - e'"* + - e''‘)(-l + A)]^! - [1 + - e^*)(-l + A - A) + A - e^*(l + A)]S' 2 ) 

B4(t) = Mgi(0)(5i - (l + A)52)g2(yo + e^**[-j/o + A^o, 


yo, 


APPENDIX B: LIST OF VARIABLES AND PARAMETERS 

Zo,i’ oxygen abundance of galaxy i\ 

Mgi'. mass of gas in galaxy i\ 
ipi’. star formation rate fo galaxy i; 

Si: star formation efficiency of gaiaxy i\ 
gas fraction of galaxy i; 
yo ■■ oxygen yield; 

R: returned fraction; 

A: outflow parameter; 

A: infall parameter; 

e: fraction of the outflowing gas from galaxy2 to galaxyl; 

Lo: outflow parameter connected to galactic fountains; 

Q.: infall parameter connected to galactic fountains; 
a: differential wind parameter. 
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